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Abstract 

We reanalyse the limits on the gravitino mass m^i2 in superhght gravitino scenarios 
derived from arguments on energy-loss during gravitational collapse. We conclude 
that the mass range IQ'^eV < 7713/2 < 2.3 x 10~^eV is excluded by SN1987A 
data. In terms of the scale of supersymmetry breaking A, the range 70GeV < A < 
SOOGeV is not allowed. 



In a wide class of supergravity models with SUSY breaking A in the TeV range, the 
gravitino can be very light: 

m3/2 = 2.5 X lQ-^eV{K/lTeVf. (1) 

Indeed, models where gauge interactions mediate the breakdown of supersymmetry , models 
where an anomalous U(l) gauge symmetry induces SUSY breaking p|, and no-scale models 
are all examples of models where a superlight gravitino can be accommodated 0. In all of 
them, the gravitino is the LSP and, furthermore, its couplings to matter and radiation are 
inversely proportional to its mass. Therefore, one may expect interesting phenomenology [^]. 
Bounds on the gravitino mass, or equivalently on the scale A have been given in the context of 
those models by various authors and have been extracted from different physical systems. In 
fact, the limits come from as distinct areas as the anomalous magnetic moment of the muon 
primordial nucleosynthesis or stellar energy drain arguments 0. In recent papers 
it has been noted that the amplitudes for gravitino processes that were used in deriving 
some of the constraints had an incorrect energy behaviour. In particular, the supernova (SN) 
bounds deduced in ref. using the effective couplings given explicitly by Gherghetta [|l^ are 
invalid as pointed out by Luty and Ponton These authors, however, when reexamining 
the limits coming from the SN1987A explosion, use an incorrect abundance of positrons in 
the core, do not discuss gravitino bremsstrahlung, and misidentify the main source of opacity 
in the stellar core. The purpose of the present paper is thus to redo the analysis that renders 
the bounds on 7713/2 (or A) following from SN collapse. Since SN considerations gave the best 
limits on A up to now ||Tl|, this is not an empty exercise. 

The relevant piece in the effective lagrangean is the derivative coupling of the goldstino x 
to photons: 

5L,ff = {e/2){M/K^fdVxF,u + h.c. (2) 

with F^j,, the electromagnetic field-strength and M is a mass that depends on the supersym- 
metry breaking model. In gauge-mediated models, M ~ mi/ Air, where mi is the left-handed 
slepton mass. Given that gravitino pairs are mainly produced via one-photon interactions, the 
sources of gravitino luminosity in stars are, in principle, gravitino bremsstrahlung in neutron- 
proton scattering, pair production in electron-positron annihilation and plasmon decay into 
gravitinos. The energy-loss rate (per unit volume) via pn pnGG is, 

(3) 
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where {p^,p)i are the 4-momenta of the initial and final state nucleons, {k^,k)i^2 are the 4- 
momenta of the gravitinos and /i^2 are the Fermi-Dirac distribution functions for the initial 
proton and neutron and (1 — fs^^) are the final Pauli blocking factors for the final proton and 
neutron. The squared matrix element can be factorised as follows, 

Y: I M;, r= {27r)W{M/AYN,Mi;j, (4) 

spins 

where Nfj_i, is the nuclear (OPE) tensor and Ggjg is the gravitino tensor in the matrix element 
squared. The factor N^^ is common to any bremsstrahlung process involving nucleons.lt ap- 
pears, e.g., in neutrino bremsstrahlung calculations and in axion bremsstrahlung calculations. 



and is given explicitly in ref fl^.On the other hand, is a tensor specific to gravitino 
bremsstrahlung. It reads, 

Gs^ = + - kiMg"" (5) 
The integration of Nfj_i, over the phase-space of the nucleons can be performed explicitly and 



the details can be found again in Raffelt's book |12|. When we contract the result with 
the gravitino tensor G3/2 and perform the integrals over gravitino momenta to complete the 
energy depletion rate, we are led to the following emissivity: 

QZfms = (8 1 92/3857r^/^)a^«^(M/A^)^Fen|T^V2/^5/2 

for non-degenerate and non-relativistic nucleons (a^ is the pionic fine-structure constant, 
is the number density of baryons, and Ye is the mass fraction of protons). However, nucleons 
are moderately degenerate in the SN core. The emissivity in the (extreme) degenerate case 
is calculated to be, 

Qf,_, = m7ryA725)a'aliM/AYpFT' (7) 

with pp, the Fermi momentum of the nucleons. Numerically, for the actual conditions of the 
star, both emissivities differ by less than an order of magnitude (about a factor of three). 
Since the actual emissivity interpolates between these two values, we shall adopt the smallest 
of the two (i.e. Q^.^^) to make our (conservative) estimates. We turn next to the annihilation 
process. 

The energy loss for the process e^{pi) + e~{p2) — > G{ki) + G{k2) can be calculated along 
similar lines as above. The spin averaged matrix element squared is in this case, 

I Mf, p= (27r)V(M/A2)%,(pi,p2)G^;2(^i,^2) (8) 

spins 

where E^y{pi,p2) equals formally the tensor 6*3^2 in eq.(5) with /ci, k2 replaced by pi,P2- The 
luminosity then is found to be, 

Q^rrn = 8a\M / A^T'e-^^^ fi'b{fi/T) /15^ (9) 
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with b{y) = {5/6)eyy-^{F^F^ + F^F^) where F^{y) = clxx"""^ / [1 + e^^*') {fi is the 
chemical potential of the electrons). The function b{y) — > 1 in the degenerate limit. Finally, 
our estimate of the plasmon decay luminosity is, 

Qp = 16C(3)a^TV(M/A2)78l7r5 (10) 

(where only transverse plasmons have been taken into account). 

Taken at face value, the bremsstrahlung rate is the largest of the three. However, Qbrems 
is overestimated since we did not consider multiple scattering effects which are present in 
a dense medium Indeed, as for the axion case[0,the gravitino bremsstrahlung rate 
probably saturates around 20% nuclear density and this should be taken into account when 
evaluating eq.(6). If we use now the values T = 50MeV, = 300MeV, and Fg = 0.3, eqs. 
(6) (with ub ^ 0.2nnuc), (9) and (10) give 

Qann : Qbrems ■ Qp ~ 1-2 X 10^ : 3 X 10^ : 1 (11) 

Therefore, a limit on A will follow from the requirement that L3/2 ~ VQann {V is the 
volume of the stellar core) should not exceed s. This constraint on the gravitino 

luminosity L3/2 imphes, in turn, 

A > 300GeV{M/A3GeVy^^{T/50MeVy^/^^{RjlOKmf^^ (12) 

or, using eq.(0), 

m3/2 > 2.3 X lO^^eV. (13) 

Of course, the previous calculation makes sense only if gravitinos, once produced, stream 
freely out of the star without rescattering. That they actually do so, for A > SOOGeV, can 
be easily checked by considering their mean-free-path in the core. The main source of opacity 
for gravitinos is the elastic scattering off the Coulomb field of the protons: 

A = 1/an = {Al^a^)Y~^p~^m^\A^/Mf (14) 

The thermally averaged cross-section for elastic gravitino scattering on electrons is roughly 
a factor Tp/rn^ smaller than that on protons and thus it does not contribute appreciably to 
the opacity. Putting numbers in eq.(14) we find: 

A ~ 1.4 X 10W(43Ge\//M)'^(A/300Ge1/)*^ (15) 

On the other hand, the calculation of Q breaks down for A < lOi^'m, i.e. for A < 220GeV , 
when gravitinos are trapped in the SN core. In this case, gravitinos diffuse out of the dense 
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stellar interior and are thermally radiated from a gravitino-sphere -R3/2- Because in this 
instance the luminosity is proportional to T^, only for a sufficiently large -R3/2 (where the 
temperature is correspondingly lower), the emitted power will fall again below the nominal 
IQ^'^erg/ s. Consequently, gravitino emission will be energetically possible, if A is small enough. 
The gravitino-sphere radius can be computed from the requirement that the optical depth 



00 



r= / drIXir) (16) 
Jr 

be equal to 2/3 aX R = Here, A(r) is given in eq.(14) with the density profile ansatz: 

p(r) = p^iRjrr (17) 
with pc = 8 X IQ^'^g/cm?, Rc = lOKm and m = 5 — 7 and which satisfactorily parameterises 



the basic properties of SN1987A ||T3|. An explicit calculation renders: 

Rs/2 = Rc[i8Yj3na^)iAyMy{m - I) / p^R^m^f^-^ (18) 
Stefan-Boltzmann's law implies for the ratio of gravitino to neutrino luminosities, 

Lg/a/L. = {R^/2/Ruf[T{R^/2)lT{R,t (19) 

where Ry is the radius of the neutrinosphere. To proceed further we use the temperature 
profile: 

T = T,{Rjrrl^ (20) 
which is a consequence of eq.(17) and the assumption of local thermal equilibrium. Now, 



taking m = 7[0, we obtain 

= {RjR,)^^l%lQYj'na'){K'/MY/p,R,m,Y'/' (21) 

By demanding that L3/2 < O.lLy and using Ry ~ 30i^m, we get 

A < lOGeV. (22) 

This in turn implies m3/2 < IQ^^eV . Since, on the other hand, the anomalous magnetic 
moment of the muon already requires ^3/2 to be larger than ~ 10"^eV^ |^, 15], we are forced 
to conclude that 

A > aOOGeV" (23) 

or, equivalently, 

m3/2 > 2.3 X IQ-^eV. (24) 
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In conclusion, we have carefully rederived the bounds on the superlight gravitino mass 
(i.e. the SUSY scale A) that follow from SN physics. These limits are completely general 
in the sense that they do not rely on other particles in a given particular model being light. 
Should other particles such as the scalar partners of the goldstino also be light, then the 
resulting bounds are necessarily tighter. In such clearly less general frame, constraints have 



also been derived in the literature [|T^] that are not subject to the criticisms mentioned in the 
beginning of this paper. They are much stronger then the ones given here and typically give 
A > 300Tel^ (or,m3/2 > 50eV ) from stellar (e.g. the Sun) evolution arguments, provided 
m3/2 < IKeV (e.g. Tq). 

Work partially supported by the CICYT Research Projects AEN95-0815 and AEN95-0882 
and the Theoretical Astroparticle Network under the EEC Contract No. CHRX-CT93-0120 
(Direction Generale 12 COMA). 
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